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The concept of thermal stability is utilized to derive expressions for the radio refractivity 
of an air parcel undergoing adiabatic compression or expansion. These expressions are of 
exponential form with scale heights of 12.5 km for the dry adiabatic process and 7.0 km for 
the wet adiabatic process. The adiabatic curves for N arc determined solely from conditions 
at the earth's surface. A new turbulence parameter, n, is derived as the difference of the 
environmental refractive index structure and the adiabatic curves for an air parcel lifted 
from the surface to ^ km. This parameter yields correlations with three hourly median 
field strength data that are not statistically different from those obtained with the much 
simpler parameter, AN. The correlations obtained with the equivalent gradient are not 
statistically different from zero. 



1. Introduction 

We shall here attempt to unify past work utilizing 
the surface value of the refractive index and the 
concept of atmospheric stability. To do so we shall 
first review the concepts of atmospheric stability 
from the viewpoint of temperature structure. In 
the process we shall show that past efforts have 
neglected the important role of the conditions at 
the earth's surface, and then we shall extend this 
work to derive an analogous expression for the radio 
refractive index. 

2. Concept of Thermal Stability 

The stratification of the earth's atmosphere by 
and large reflects the control of the earth's gravita- 
tional field. This average structure is, however, 
constantly disturbed by the intrusion of thermal 
plumes of heated air rising from the differentially 
heated surface of the earth as well as the mechanical 
mixing produced by the passage of air currents. It 
is customary to assume that parcels of air are forced 
upwards through the normal stratified atmosphere 
without mixing with the environmental air mass. 
Since the parcel is forced to rise, be it by convection 
or mechanical mixing, it is also assumed that it 
expands and cools without exchange of heat; i.e., 
adiabatically. In such a process it can be shown 
that the pressure, P, is given by 
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where T is the temperature (°K), g, the acceleration 
of gravity; R } the gas constant for air; and a, the 
constant lapse of temperature with height is 



-a=+9.8°C/km. 
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The zero subscripts in (1) refer to conditions at the 
earth's surface. Equations (1) and (2) are taken to 
hold as long as saturation does not occur. It is 
further assumed that the mixing ratio 
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is constant throughout the process. In (3) e is the 
partial pressure of water vapor. 

A parcel of air forced to rise will cool to the dew- 
point if forced aloft far enough. Further rising will 
lead to condensation. Conversely a parcel of air 
which sinks down will have work done upon it adi- 
abatically and will become warmer than it was in 
its original elevated position. Air following a con- 
densation curve has a value of a, a*, of 6 °C/km. 
In addition, as water vapor condenses out of the 
parcel, the vapor pressure will decrease. We here 
assume that the vapor pressure will follow that of 
the saturation vapor pressure curve, e s , as given by 



e sz = e sl exp [— 0a* (z — I)], z>l, 



(4) 



where I is the height of the lift condensation level 
and /3=0.064( o C)-\ The value of I is given [Hewson 
and Longley, 1944] with sufficient accuracy for prac- 
tical applications by 
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where the zero subscript indicates that the difference 
between the temperature and the dewpoint, T d , need 
be evaluated only at the earth's surface. 

Consider now a parcel rising through some envi- 
ronmental distribution of temperature. The dry 
adiabatic lapse of temperature with height is greater 
than that observed on the average in the atmosphere 
(6 °C/km) and thus the parcel becomes cooler, and 
more dense, than the environmental air and will sink 
to its initial conditions with the removal of the lift- 
ing force. If, on the contrary, anywhere in its tra- 
jectory it becomes warmer, and thus less dense, it 
will become unstable and rise of its own buoyancy 
through the environmental air. Past radiometeor- 
ological studies have taken the area between the adi- 
abatic curves and the environmental temperature 
distribution on a pressure (~ height) versus tem- 
perature chart as a measure of atmospheric stability. 
On the general argument that the effects of atmos- 
pheric stability, or turbulence, influence radio waves 
only through the radio refractive index, we shall ex- 
tend the above concepts to the radio refractive index. 

3. Adiabatic Lapses of N 

The radio refractivity, with the usual units, is 
given by 
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Introducing the dry adiabatic variations of P, T, 
and e as given above one obtains the dry adiabatic 
decrease, N d , of N: 

which is valid for z<l. The distribution of N d with 
z may be evaluated by taking the logarithms of both 
sides of (7), 

In *W. 2 V „ + (X- 2 ) 1 n( 1 _|) +h (l-^) 

<S) 
and expanding the logarithms on the right-hand side 
(r.h.s.) of (8) with the neglection of second order terms 
where appropriate for our applications of z<3 km: 
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When evaluated for the standard conditions of 
T =288 °K (15 °C) and P = 1013.2 mb, T d becomes 



r d =0.08 (km)" 



(12) 



It is sometimes more convenient to use the dry 
adiabatic scale height, H d ; 



ff d =^=12.5km. 



(13) 



Above the local condensation level, the expression 
for N becomes that of the wet adiabatic, N w : 

+W l exp ((3a*z) V (14) 

where, for convenience 

A =77.6£andl^=™(f^. 
J- 1 J- 1 

(Throughout this development the subscript I refers 
to the local condensation level.) Applying the same 
approximations as in the derivation of (9) one obtains 
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For the standard conditions assumed above and e 
corresponding to 60 percent relative humidity H w 
~7.5 km. Again, the form of the wet adiabatic 
lapse of N is given by 



N w =N t exp (—T w z), z>l. 



(17) 



We shall now apply these results to the derivation 
of a refractive index turbulence parameter. 

4. Turbulence Parameter, n 

Analogously to the concept of thermal stability, 
we define II as the area between the environmental 
N(z) curve and the appropriate adiabatic decrease 
of N. That is 
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where the integration is arbitrarily taken from the 
surface to 3 km. This then extends the integration 
well above the crossover height of the radio horizon 
rays from transmitter and receiver normally encoun- 
tered in tropospheric propagation. The integral 
for II may be written 

n= CncIz- f V exp (-T d z)dz 
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N t exp [-T a (z-l)]dz. (19) 
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Noting that 

N t mN a exp(-Tj), 

(19) then becomes, upon integration, 

n=fV(/z+^[exp(-lV)-l| 

JO -L d 
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(20) 



II is thus'depeadent upon both N profile character- 
istics and initial conditions. In fact, all the terms 
on the r.h.s. of (20) save 



/; 



Ndz 



are determined from the conditions at the earth's 
surface. One might expect, then, that II would 
constitute a useful radiometeorological parameter 
since it includes the parameter JV" , the concept of 
atmospheric stability, and the integrated iV-profile 
characteristics. We shall now apply this new 
parameter to the radio data presented in detail in 
part II. The values of H d and H w utilized are 
those for standard sea-level conditions; i.e., 12.5 km 
and 7.5 km respectively. 

5. Comparison of Radiometeorological 
Parameters 

The various radiometeorological parameters dis- 
cussed in part I plus II as derived from the Rantoul 
and Joliet radiosondes were tested for their relative 
utility by comparison with the radio data of WBKB- 
TV and WENR-TV (see table 2, part II). Correla- 
tion coefficients between the radio data (three hourly 
medians) of each path and the various radiometeoro- 
logical parameters were determined for each of the 
profile categories given in part II, table 3, as well as 
for all available data. It was not possible to fit all 



of the data into the various categories since these 
were deliberately made very strict in order to explore 
the utility of each parameter within each of the 
propagation mechanisms assumed for the profile 
categories (pt. II, table 1). Thus it is observed that 
the sum of the individual categories does not equal 
all available data. The results of these correlations 
are given in table 1. Before discussing these results 
it is well to note that the crossover height of the 
tangent rays from each end of the path were cal- 
culated by detailed ray tracing [Bean and Thayer, 
1959] for each radiosonde ascent. The values as 
determined from the Rantoul and Joliet soundings 
were averaged. This average crossover height was 
then used to determine the values of II below the 
crossover height (II"), above the crossover height 
(IT), and g e . As listed in table 1, II represents the 
value of (18) for the height increment zero to 3 km. 
The conclusions reached below are based only upon 
the data, of category six (all data) of table 1; the 
sample size of the individual categories sometimes 
being so small as to prohibit the attachment of any 
significance to the observed results. It is thus evi- 
dent from table 1 that the results obtained with II 
are comparable with those using II', due supposedly 
to the crossover heights being on the order of 100 m 
for these paths. II" appears to contribute little to 
the final value of II and is poorly correlated with the 
radio data. The overall conclusion that one reaches 
from the data of table 1 is that the most promising 
of the parameters considered is II closely followed 
by AN and N s . In fact, the difference between the 
results obtained with II and AN gives one pause to 
consider if the added complexity of determining II 
is worth the relatively small gain in correlation with 
the radio data. Further, examination of the 95 
percent confidence belts of the correlations [Crow, 
Davis, and Maxfield, 1960] shows that the results 
obtained with AN and II are not statistically different 
within this confidence band. Similarly, it is concluded 
that g e yields correlations not statistically different 
from zero. 



Table 1. Summary of correlations 















Radio meteorological variable 




Station 


Category* 


Sample 
size 


































n 


n' 


n" 


0e 


Qo 


Ns 


AN 


WENR-TV (179.75 


1 


68 


-0. 379 


-0. 380 


-0. 186 


0.189 


-0. 188 


0.141 


-0. 346 


Mc/s) 


2 
3 


21 

101 


-.206 
-.175 


-.268 
-.175 


.200 
.009 


.004 
.029 


.135 
-.030 


.053 
-.094 


-.265 




-.455 




4 


6 


-.363 


-.350 


-.432 


.499 


-.415 


.090 


.087 




5 


196 


-.566 


-.573 


-.031 


.056 


-.177 


.387 


-.564 




6 


518 


-.462 


-.460 


-.017 


.072 


-.176 


.289 


-.454 


WBKB-TV (71.75 




















Mc/s) 


1 
2 


67 
27 


-0. 295 
-.207 


-0. 298 
-.184 


-0. 223 
-.138 


-0. 093 
-.104 


0.079 
.094 


0.011 
. 168 


-0. 169 




-.161 




3 


113 


-.269 


-.271 


.076 


.165 


-.081 


.070 


-.349 




4 


6 


-.244 


-.220 


-.685 


.308 


-.369 


.257 


.348 




5 


213 


-.552 


-.556 


-.058 


.062 


-. 115 


.417 


-.464 




6 


564 


-.438 


-.440 


.009 


.043 


-.103 


.213 


-.350 



*Category 

l = Unstra titled. 

2= Ground based duct. 

3 = Elevated layer. 

4=Subrefractive layer. 

5=1 through 4 combined. 

6= All data for period of record. 
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The above data also allows one to evaluate a 
parameter of the form 

H l =ag e +bU r 

which is analogous to Misme's recently suggested 
M parameter [Misme, I960]. The parameter Hi 
given above is physically more realistic for radio 
purposes since it incorporates a measure of the 
refractive index turbulence rather than thermal 
stability. The correlations obtained between field 
strength and Hi are 

WENR-TV: 0.46 
WBKB-TV: 0.44 

which represents, for the significant figures carried, 
no improvement over the use of II alone. This last 
analysis suggests that it might be desirable to con- 
sider other combinations such as 

H 2 =a 2 AN+b 2 n. 

The correlations for this case are 

WENR-TV: 0.49 
WBKB-TV: 0.44 

which represents a slight improvement over II alone. 



6. Conclusions 

In conclusion, an atmospheric stability parameter, 
II, based upon the radio refractive index has been 
developed. The IT parameter, like thermal stability, 
is an integrated effect and thus many combinations 
of vertical temperature and humidity distributions 
may yield the same value of II (the same is true of 
thermal stability). Nevertheless, the correlation 
obtained between II and radio data on a three hourly 
median basis is perhaps the most encouraging ob- 
tained to date but does not yield correlations signifi- 
cantly different from those obtained with the much 
simpler parameter, AN. This conclusion could well 
change with broader experience on other paths and 
radio frequencies. 
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